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Upward  Flame  Spread  on  Vertical  Surfaces 


1.0  INTRODUCTION 

In  this  report,  a  vertical  flame  spread  model  is  developed  and  validated  against  literature  data.  This 
work  is  a  part  of  the  U.S,  Navy  Materials  Test  program  and  is  designed  to  provide  a  technique  for  specifying 
the  performance  required  for  tte  U.S.  Navy  material  applications  in  terms  of  small-scale  tests.  This  program 
will  result  in  known  fire  performance  of  material  applications  and  will  allow  manufacturers/developers  to 
provide  cost-effective  materials  with  the  required  performance,  leading  to  cost  saving  and  known  performance. 

The  upward  flame  ^read  model  on  vertical  surfaces  is  part  of  the  modeling  effort  designed  to  develop 
a  general  modeling  framework  i^iiich  will  assess  performance  of  materials  in  Navy  fire  scenarios  firom  small- 
scale  test  data.  This  model  is  formulated  on  the  basis  of  a  review  of  the  fire  dynamics  literature  relevant  to 
fire  growth  presented  by  Williams  and  Beyler,  1994  (reference  (1)). 


2.0  THEORY  AND  DESCRIPTION  OF  THE  MODEL 

A  con^uter  model  has  been  developed  which  calculates  the  flame  spread  on  a  vertical  wall  subjected 
to  a  line  igmtion  source.  Flame  spread  is  calculated  using  sub-models  derived  from  the  literature  utilizing 
u^Hits  determined  from  cone  calorimeter  tests  (reference  (2)).  The  model  is  formulated  based  on  existing  sub¬ 
models  for  (1)  ignition,  (2)  material  heating,  pyrolysis,  and  burning  rate,  (3)  flame  spread,  and  (4)  flame  and 
surface  heat  transfer.  The  details  for  each  component  of  the  analysis  will  be  described  in  the  following 
sections. 

2.1  Model  Summary 

The  computer  model  calculates  the  flame  spread  on  a  vertical  surface  by  breaking  up  the  surfoce  into 
a  large  number  of  elements.  The  conditions  of  each  element  are  independently  computed.  The  centroid  of 
an  element  is  assumed  to  be  representative  of  the  entire  element  (Fig.  1).  There  are  a  number  of  global 
conditions  that  are  calculated  by  summing  the  contributions  from  each  element:  the  heat  release  rate,  the 
height  of  the  pyrolysis  region,  the  flame  height  above  the  base  of  the  wall,  and  the  height  of  burnout  front 
above  the  base  of  the  wall. 

Each  element  is  in  one  of  four  states:  (1)  preheat  (above  the  flame),  (2)  preheat  (exposed  to  the  flame), 
(3)  burning,  and  (4)  consumed.  The  model  keeps  track  of  these  conditions  for  each  element,  and  the  model 
stops  when  either  the  user  entered  simulation  time  is  reached,  the  fuel  is  entirely  consumed,  or  the  flampi 
propagation  ceases. 
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Figure  1  Line  fire  configuratiQn  against  the  wall 


Wall  Flame  Heights  and  Heat  Fluxes  to  the  Surface 

In  the  present  model,  the  wall  can  be  heated  by  three  sources:  (1)  an  imposed  external  heat  flux,  (2) 
a  line  fire  pla^  against  the  base  of  the  wall  (the  ignitor),  and  (3)  the  wall  flame  itself.  It  has  been  shown 

that  the  heat  flux  to  the  waU  from  line  fires  against  the  wall  and  the  wall  itself  can  be  correlated  in  the  same 
manner. 


^*8-  2  shows  a  correlation  of  the  heat  flux  from  a  methane  line  burner  to  an  adjacent  wall  along  with 

similar  results  for  fires  involving  liquid  fuels  saturated  on  the  lower  portion  of  a  wall.  This  correlation  can 
be  expressed  as  follows; 


20—  ^<0.34 


=  6.23y‘‘‘”  — 

0.34  s  ^  <  0.7 

yf 

=  3.59  y  — 

0.7 

V. 

(1) 


where  q"  is  the  incident  heat  flux  from  the  line  burner  to  the  wall  surface  (kW/m^), 

y  is  the  height  above  the  base  of  the  wall  (m),  and 
yf  is  the  flame  height  above  the  base  of  the  wall  (m). 

This  correlation  is  usefiil  for  representing  methane  burner  ignitor  and  commonly  used  liquid  ignitor  fuels. 

An  alternative  approach  available  in  the  model  is  the  use  of  one  of  the  waU  fire  heat  flux  correlations 
fisted  later  m  this  section.  These  correlations  are  based  on  more  luminous  fuels  which  are  typical  of  waU 
flames  and  non-methane  line  burners.  The  result  of  using  the  wall  heat  flux  correlations  is  a  more  rapid 
igmtion;  however,  the  fire  growth  will  be  similar. 

Equation  (1)  requires  the  determination  of  the  flame  height.  There  are  several  flame  height 
TOrrelations  available  for  a  line  fire  against  a  wall  or  wall  fires  (references  (4)  and  (5)).  Each  of  these  takes 
the  form: 


Yf  -  C-  0^^  for  —  >  R 

y 

2^ 

C 

yf  -  R 

y 


(2) 


\^ere  yj  is  the  wall  flame  height  (m), 
yp  is  the  pyrolysis  height  (m), 

Q'  is  the  heat  release  rate  per  unit  length  of  the  line  burner  (kW/m), 
Cf  is  a  constant  (m®'’/kW“),  and 

R  is  the  minimum  ratio  of  the  flame  height  to  the  pyrolysis  height. 
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Fig.  2  -  Wall  hMt  flux  distribution  due  to  wall  fires.  The  solid  symbols  are  for  a  lina  burner 

fire  adjacent  to  a  wall  and  the  open  symbols  are  for  liquid  saturated  wall  fires  with  low  flatw 
intensity  (reference  (3)) 
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Delichatsios  (reference  (4))  recommends  Cr  =  0  052  anrf  Tu  nnH  /  <• _ 

TOommoKl  C,  -  0.0666  m»/kW”.  The  eet<nd  pm  of  Eq J^)  S^Jiceble  »  weU  to 

height  to  ^^®*8ht.  Based  on  work  by  Tewarson  (reference  (6)),  R,  the  mininnim  ratio  of  the 

height  to  the  pyrolysis  height,  is  in  the  range  of  l.O-l.l.  One  of  the  goals  of  the  comparison  with 

»hich  velue  of  C,  fiom  to  Utotoe  i,  most  sppmprim  for  us.  m  vetod 

‘“^°sity  of  the  flame  may  increase  over  that  typical  of  the 

^  ^  “se  of  a  different  heat  flmfcorrelattoS 

^r  lotion  of  the  w^.  The  mcident  heat  flux  from  wall  flames  to  the  surface  has  been  exuerimentallv 
etermmed  for  a  numbw  of  materials.  Fig.  3  shows  wall  heat  flux  distributions  from  wall  flameMi  to  the 

“sf ‘t?  «"dos«  u^t  of  to  stody  stete  to  flux  to  ^sto  as^ZSS 
3a  ^  been  developed  as  shown  in  Fie. 

flames  back  to  toe  Mitler  (reference  (8))  also  calculated  heat  fluxes  from  toe  wall 

An  upper  bound  fit  to  Fig.  3, 


=  12.7 


(y-yj) 

(yf-y^) 


(y-y») 

(yf-y^) 


9".  =  30 


(y-yj) 

(yf-yp 


a  lower  bound  fit  to  Fig.  3, 


=  2.71 

*  ✓  V 


Oy-n) 


(y-yj) 

(yry„) 


k  0.43 


J//  .  to  kW 

1  «  -  1*  —7 
m 


(y-y>) 


<  0.43 
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Wall  Flame  Heat  Flux  (kW/irf) 


Figure  3  Rame  heat  flux  distribution  for  wall  fires  in  terms  of  y/yj . 
(reference  (7)) 


100 


10 


0.1 


An  if^por  bound  fit  Fi9«  3: 


fitFiivtJ: 


«*.  -  2J  . 
a  toMT  beiBd  fit  Fipid  3: 


tyy>> 


<0.7 


*2::ZtL,o 

<y«y>) 

<y/-y*j 

<yy>)  , 

tyr-y*i 


\V'^ 


0.01 
04101 


0.01  Oil  1 

yhft 

Rgure  3a  Wall  flame  heat  flux  carrelations 
(reference  (8)) 
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and  an  average  fit  to  Fig.  3; 


=5.39 

®  .Tv 


(y^-y,) 


(y-yt) 

(yf-yt) 


^0.53 


=  23  ^ 


Cy-n) 

Oy-y*) 


<0.53 


where  y  is  the  height  above  the  base  of  the  burning  wall  (m), 

Yb  is  the  height  of  the  burnout  front  (m),  and 

yr  is  the  height  of  the  flame  tip  above  the  base  of  the  burning  wall. 

Under  spreading  conditions  with  burnout  occurring,  the  yj  should  be  replaced  by  yj  -y^, 

The  selection  among  these  correlations  wUl  depend  on  the  flame  luminosity  and  the  level  of  conservatism 
required.  As  tht  various  heat  flux  correlations  available  span  a  considerable  range  of  incident  heat  fluxes 
one  of  the  goals  of  the  comparison  with  experimental  data  is  to  determine  which  of  these  correlations  from 
the  literature  is  most  ^propriate  for  use  in  vertical  flame  spread  predictions. 

2.3  Surface  Heating  and  Ignition 

Prior  to  ignition,  the  waU  surface  is  assumed  to  be  a  semi-infinite  one-dimensional  slab  with  a  time 
di^en^nt  surface  heat  flux  determined  from  the  waU  flame  and  external  sources  as  described  in  Section 
conduction  model  us^  in  this  computer  model  is  an  approximate  solution  to  the  semi-infinite  slab 

temperature  profile  and  an  integral  solution  (reference  (9)).  TOs  method 
was  selected  based  on  its  excellent  predictive  performance  and  computational  efficiency. 

(refer  ^  *>e  approximated  by  the  following  equation 

Tiz)  =  +Bz^ +Cz +D  (9) 

where  z  is  the  distance  from  the  surface  and 

the  constants  A,  B,  C,  and  D  may  be  determined  from  the  boundary  conditions: 


at  z  =  6,  T  =  T 


at  z  =6, 
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at  z  =  0 


(9d) 


The  resulting  temperature  profile  becomes 


4"  6 
T  =T  +1-SS— 

*  "  Zk 


1-- 


(10) 


>^ere  T,  is  the  initial  surface  temperature  ("C  or  K), 

6  is  the  thennal  penetration  depth  (m), 
z  is  the  depth  into  the  wall  surface  (m),  and 
k  is  the  thermal  conductivity  (kW/m-K). 

This  profile  may  be  inserted  into  an  integral  form  of  the  energy  equation: 

-.(I) 


(11) 


M^ere  a  is  the  thermal  diflusivity  (mVsec). 
differential  equation: 


The  surface  temperature  at  any  time  is  then  given  by  the  foUowing 


d_ 

dt 


( -  K  Y 


mt 


3  kgs 


(12) 


v^cre  p  is  the  density  (kg/m^)  and  c  is  the  thermal  heat  capacity  (kJ/kg-K), 
of  the  wall  surface  temperature  is: 


The  net  heat  flux  as  a  function 


-  oei7;(o^-?^j  -hjj'j.i)  -  rj 

where  q\  is  the  incident  heat  flux  firan  the  flame  to  the  surface  (kW/m*), 
q"^  is  the  externally  imposed  radiant  heat  flux  (kW/m^, 

h,  is  the  heat  transfer  coefficient  (kW/m*-K)), 

T,  is  the  surface  temperature  (®C  or  K),  and  * 

T,  is  the  ambient  temperature  (“C  or  IQ. 


T  <*^8  preliMting  prior  to  flame  initiation  where  flamefphime 

of  Equation  (12)  and  diffetentiating  EqLon 
(13)  with  respect  to  time  yields  a  differential  equation  with  T,(t)  as  the  only  unknown  term: 
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(14) 


2 


dr, it) 

dt 


q"jt)  dt  3  ^pc 


dt 


=  -  4  T,{tf 


dt 


(15) 


Equations  (14)  and  (15)  arc  solved  in  the  flame  spread  model  using  a  fourth  order  Runge*Kutta  method  for 
each  time  step.  This  provides  a  very  efficient  method  for  predicting  the  heating  of  the  surface  up  to  the 
ignition  temperature.  The  thermal  inertia,  kpc,  is  determined  from  ignition  experiments  in  the  cone 
calorimeter. 

In  the  early  stages  of  the  temperature  rise,  the  differential  equations  do  not  converge  due  to  a 
singularity  i^en  T,  (t)  is  equal  to  T^.  For  small  temperature  rises  Equation  (16),  a  lower  order  method  without 
this  singularity,  was  employed  to  determine  the  temperature  rise  of  the  surface.  This  method  uses  an  average 
net  surface  heat  flux  and  is  a  single  calculation  from  the  initial  temperature: 


N  « 


(16) 


where  q"„^  is  the  average  net  heat  flux  received  at  the  surface  (kW/m  ^ )  over  time  t  (sec).  It  was  determined 
that  Equations  (14)  and  (15)  have  convergence  problems  up  to  a  surface-initial  temperature  rise  no  greater 
than  2-3  K  and  that  both  methods  wctc  in  fair  agreement  up  to  a  temperature  rise  of  50-100  K.  Equation  (16) 
is  used  in  the  model  for  only  one  or  two  time  steps.  The  material  properties,  the  heat  flux  level,  and  the 
duration  of  the  heat  flux  levels  all  determine  where  the  integral  method  begins  to  converge. 

2.4  Heat  Release  Rate 

Once  an  element  has  reached  the  ignition  temperature,  it  is  allowed  to  pyrolyze,  bum  and  contribute 
energy  to  the  wall  fire.  The  heat  release  rate  of  the  wall  is  determined  by  summing  up  the  heat  release  rates 
of  each  burning  element: 


n 


=  m"  (x,  y,  z,  t)  A  {x,y) 

/•I 


(17) 


wfrere  n  is  the  number  of  elements  ^^ch  are  burning  at  time  t  (sec), 

AHe  is  the  heat  of  combustion  of  the  wall  material  (kJ/kg), 

ih"(x,  y,  z,  t)  is  the  spatially  and  temporally  varying  mass  loss  rate  per  unit  area  (kg/m^),  and 
A(x,  y)  is  the  area  of  the  element  under  consideration  (m^. 
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Both  elemeotal  area  and  the  mass  loss  rate  are  allowed  to  vary  from  element  to  element,  althn^igh 
mere  is  curren^  no  convenient  way  to  specify  elements  of  variable  size  implemented  in  the  current  model, 

i.e.,  adaptive  gridding.  The  mass  loss  rate  for  each  element  is  given  by 


rn  "iO  = 


(18) 


A^ere  Ahg  is  the  heat  of  gasification  (kJ/kg)  as  determined  from  the  cone  calorimeter.  The  net  energy  (q"n-) 
absorbed  by  each  burning  element  can  be  calculated  from  the  following  equation; 


4'Lp  aeir^iO*  -  r*]  ,  (19) 

Q"in>  is  the  radiant  imposed  heat  flux, 

T,  is  the  surface  pyrolysis  temperature  (°C  or  K), 

To  is  the  ambient  temperature  (°C  or  K), 

o  is  the  Stefan-Boltzman  constant  (5.668E-08  W/m^-K^,  and 

€  is  the  surface  emissivity. 

2,5  Material  Property  Data 

The  material  properties  required  by  the  model  are  determined  using  the  cone  ralnrimyiter  ASTM 
standard  test  method  E-1354  (reference  (2)).  These  properties  are  as  follows: 

1.  Ignition  temperature,  Tig,(K); 

2.  Thermal  inertia,  kpc,  (kW/m^-K)^  sec; 

3.  Effective  heat  of  gasification,  Ahp  (kJ/kg);  and 

4.  Heat  of  combustion,  AH^,  (kJ/kg). 

The  ignition  temperatoe  m^  be  measured  directly  in  the  cone  calorimeter  test  or  it  mjy  be 
estimated.  If  the  value  is  estimated,  a  set  of  ignition  temperature-thermal  inertia  values  will  result  and 
additional  information  is  necessary  in  order  to  determine  the  most  appropriate  pair.  The  thermal  inertia  is 
determined  by  trial  and  error  as  follows: 

The  ignition  times  for  a  series  of  heat  fluxes  is  first  obtained  from  the  cone  for  the  material  in  question.  The 
flame  spread  model  is  then  used  to  predict  the  igmtion  time  for  the  given  ignition  temponture,  applied  radiant 
flux  and  th^al  inertia.  By  holding  the  ignition  temperature  constant  and  altering  the  thermal  inertia  such 
4at  &e  ignition  times  ^e  predicted  for  each  applied  radiant  flux,  a  thermal  inertia  is  established  for  the 
ignition  temperature.  Figures  5  and  6  show  the  actual  and  predicted  igmtion  times  for  various  flux  levels 
at  two  assumed  ignition  temperatures,  3 17°C  and  350°C. 

The  effective  heat  of  gasification  is  derived  from  the  cone  calorimeter  mass  loss  rate  measurements 
(Fig.  7).  The  heat  of  gasification  varies  with  time,  but  an  average  value  may  be  used.  Use  of  heat  of 
gasification  that  varies  substantially  from  instant  to  instant  (as  is  often  the  case  in  cone  measured  data)  would 
produce  he^  release  rate  and  flame  height  predictions  with  oscillations  of  similar  intensity.  The  hrat  of 
gasification  is  estimated  from  the  following  equation: 
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Ignition  Time  (sec) 


350 


Incident  Heat  Flux  (kW/iii2) 

Hguie  5  Con^janscui  of  ©qperimental  and  predicted  igmdcm  time  for  plywood  at  Tj  =  317  ®C,  and 
kpc  =  0.73  (kW/n^  K)^  sec.  Data  fiom  Quintiere  and  HaiMeroad,  1985  (reference  (10) ) 


.2 


Ignition  Time  (sec) 


Measured  ignition  time 


Predicted  ignition  time 


Incident  Heat  Flux  (kW/m2) 


FigUTB  6  Comparison  of  csqjerimental  and  predicted  ignititm  timft  for  nlvwDod  at  T-  =  350  ®r'  anH 

~  0.475  (kW/m?  sec.  Data  from  Quintiere  and  Harideroad,  1985  (reference  (10) ) 


Heat  of  Gasification  Ah  (kJ/g) 


Nondimensional  Pyrolysis  Depth  (5/5^) 


Figure  7  Cone  Calorimeter  heat  of  gasification  for  cardboard  at  various  heat  flux  exposure 
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(20) 


^1(0 

rh'ii) 


where  q  IS  the  net  heat  flux  per  unit  area  incident  on  the  surface  of  the  material  and  m"  is  the  mass  loss  rate 
per  xmt  area  The  net  heat  flux  is  the  summation  of  the  imposed  radiant  flux,  the  radiative  and  convective 
tlux  the  flame  to  the  surface  and  the  radiant  heat  loss  from  the  surface  to  the  surroundings  These  terms 
are  listed  m  Equation  19  of  the  previous  section.  In  order  to  generalize  the  heat  of  gasification  values  for 
variable  heat  fl^  levels,  the  time  dependent  heat  of  gasification  is  transformed  to  a  thickness  dependent  heat 
ot  gasification  by  irplacing  time  with  the  pyrolysis  penetration  depth  at  that  time  (6): 


AA,(6)=  AA^[(5(r)l=  A/,^(0 

where  6(t)  is  the  nondimensional  pyrolysis  depth  at  time  t  and  is  determined  from 

n 

-  ntj. 


(21) 


6(/)  = 


_  <=i 


(22) 


ntf  -  ntj 


At,Js  &e  tinM  interval  of  the  mass  loss  rate  measurement  (s),  A  is  the  specimen  surface  area  (m^ 

nif  IS  the  final  specimen  mass  (kg),  and  lUi  is  the  initial  specimen  mass  (kg). 


c  1  “™*’“stion  is  calculated  in  accordance  to  the  methods  given  in  ASTM 

h  (reference  (2)).  This  is  given  as: 


Aft  -=i:i 

c,  eff 


(23) 


-  nij 


where  4  is  the  heat  release  rate  at  time  i  as  determined  by  the  test 


3.0  PROGRAM  DESCRIPTION 

.  j  program  is  written  in  LAHEY  FORTRAN  90,  and  the  foUowing  is  a 

detailed  description  of  the  main  program  and  the  major  subroutines. 


There  are  subroutines  for  reading  data,  calculating  results,  and  reporting  the  results  to  a  file  (x-  printer. 
The  mam  program  calls  the  subprograms  or  subroutines  that  read  the  data  and  the  surfiice 

temperature  at  each  element/node,  heat  release  rate,  and  flame  height  The  subroutines  structure  of  the  model 
IS  illustrated  in  the  flow  chart  of  Fig.  8. 
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3.1  Sub-models  in  the  Numerical  Simulation  Computer  Program 

The  program  structure  is  divided  into  a  main  driver  responsible  for  calling  the  various  subroutines 
in  the  required  order:  Initialization  subroutines  that  are  used  only  once,  the  calculation  and  updating 
subroutines  and  junctions,  and  the  output  and  termination  subroutines.  The  main  driver  calls  three 
initialization  subroutines  during  program  start-up:  Getinput,  DataCheck,  and  Initialize.  Initialize  will  call 
either  Burn_initial  or  Pool_initial  depending  on  the  ignition  source  exposing  the  wall.  Getinput  reads  the 
data  from  an  input  file  called  ProjecLDat.  This  file  contains  all  of  the  needed  material  properties,  the  wall 
dimensions,  analysis  flags,  and  time  and  time  step  data.  There  are  several  analysis  flags  available.  There  can 
be  a  cone  calorimeter  simulation,  the  ignition  source  can  be  entirely  separated  from  the  heat  release  rate  and 
flame  height  calculations,  several  debug  settings  exist  that  dump  a  number  of  additional  variables  to  output 
files,  and  the  ignition  source  type  can  be  selected.  Getinput  is  also  responsible  for  reading  in  the  heat  of 
gasification  data  from  the  file  IIeat_Gas.Dat.  The  subroutine  DataCheck  checks  the  input  data  to  ensure  that 
there  are  no  errors  that  could  potentially  cause  the  program  to  lock  up.  If  an  error  is  detected,  such  as  a 
negative  temperature  or  an  invalid  flame  height  correlation  selection,  the  program  is  aborted.  Next  the  driver 
calls  Initialize  subroutine.  This  subroutine  is  responsible  for  setting  the  initial  values  for  all  the  variables. 
Depending  on  the  ignition  source,  this  subroutine  will  call  either  Burnjnitial  or  Pooljnitial  to  further 
initialize  variables  specific  to  the  ignition  source. 

The  main  calculation  loop  in  the  main  driver  is  a  triple  loop  in  time  and  the  x  and  y  coordinates. 
After  the  first  time  step  the  driver  calls  the  subroutines  GetHG,  mlr,  SetFlags,  and  SetFlameHdght.  These 
subroutines  update  the  status  of  each  element  and  calculates  the  current  flame  height  on  the  wall.  Next  the 
driver  calls  the  subroutine  flux  which  determines  which  surface  temperature  calculation  subroutine  to 
Secant  or  Runge-Kutta,  depending  on  the  overall  temperature  rise  of  the  element.  After  all  the 
surface  temperatures  have  been  updated,  the  driver  calls  EndCheck  to  determine  if  either  the  fire  has 
extinguished  or  all  the  mass  has  burned  out.  If  either  is  true,  subroutine  OutPut  is  called  and  the  program 
terminal.  Otherwise  the  time  is  incremented  and  the  procedure  repeats.  The  following  is  a  listing  and  brief 
description  of  the  subroutines,  functions,  and  input  and  output  files  associated  with  the  flame  spread  program. 

3.1.1  Subroutines  and  Program  Driver 

ntainQ  -  main  driver  responsible  for  calling  subroutines  in  order  and  incrementing  time  and  spatial  loops. 

GetInputQ  -  Reads  ProjeetDat  and  heat_gas.dat  input  files  and  stores  the  data  in  common  blocks. 

DataCheck(integer  errno)  -  Checks  the  input  data  for  errors  that  may  crash  the  program  and  if  an  error  is 
detected,  returns  the  error  number. 

InitializeO  -  Initialize  variables  common  to  all  exposures. 


BurnJnitialO  -  Initialize  variables  common  to  only  a  line  burner  ignition  source. 

Pool_InitialO  •  Initialize  variables  common  to  only  a  pool  fire  ignition  source. 

SetFlags(real  hr)  -  Adjusts  the  flag  state  for  each  element.  Elements  can  be  in  preheat  above  the  flame,  in 
preheat  behind  the  flame,  ignited,  burned  out,  or  inert  (void). 

EndCheckQogical  mflag)  -  Determines  if  either  the  wall  material  has  been  completely  consumed  or  if  the 
fire  has  burned  out  and  returns  a  .TRUE,  or  .FALSE. 
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Output(real  hrr,  integer  icount)  -  Writes  the  necessary  output  data  to  flame.dat  (flame  height  and  pyrolysis 
data)  and  heatreal.dat  (heat  release  rate  data).  If  debug  mode  is  invoked  additional  files  include  roords  dat 
(elemental  data)  and  temps.dat  (elemental  temperature  data)  and  flux.dat  (elemental  flux  data). 

Secantfmteger  ni,  integer  nj,  integer  icount,  real  qa^,  real  rtime)  -  Solves  for  the  surface  temperature  at 
time  'rtime'  using  the  constant  heat  flux  equation.  This  is  used  only  for  a  specified  temperature  difference 
fi^om  the  initial  surface  temperature  at  time  0  to  avoid  numerical  instability. 

Runge-Kuttafmteger  ni,  integer  nj,  real  qoq>)  -  Solves  the  cubic  temperature  profile  differential  equation 
for  the  surface  temperature  at  time  'rtime'. 

flux(integ^  kflag,  integer  ni,  integer  nj,  real  rtime,  integer  icount)  -  Determines  the  exposure  flux  for 
element  (ni^j)  and  calls  the  respective  surface  temperature  subroutine. 

GetHG(real  heatjgas,  integer  kflag,  integer  ni,  integer  nj)  -  Updates  the  current  heat  of  gasification  value 
for  element  (ni^j)  based  on  the  pyrolysis  penetration  depth. 

mlr(real  heat_gas,  integer  kflag,  real  rmassloss,  integer  ni,  integer  nj)  -  Computes  the  mass  loss  rate  per  unit 
area  for  element  (ni,nj). 

KiUBurnerO  •  Shuts  off  the  burner  immediately. 

TurnBurnerOff(real  rtime)  -  Shuts  off  the  burner  after  all  the  material  behind  the  initial  burner  flame 
burned  out. 

SetFlameHeightfinteger  nend)  -  Adjusts  the  current  flame  height  based  on  the  updated  total  heat  release  rate 
per  unit  length.  If  there  is  no  heat  release  rate  (and  consequently  no  flame  height),  it  returns  a  flag  inrfiraring 
the  program  should  terminate. 

GetBurnHeat(real  rtime)  -  Ckrmputes  the  burner  heat  release  rate  as  a  linear  interpolation  of  the  final  burner 
heat  release  rate  over  some  specified  time  to  achieve  this  mavinnim 

3.J.2  Functions 

real  GetFlux(real  ycGF)  -  Calculates  the  heat  flux  for  an  element  at  a  height  ycGF  using  one  of  several  heat 
flux  correlations. 

real  q_bum(ycQB)  •  Calculates  the  burner  contribution  to  the  surface  incident  heat  flux  if  the  burner  fire  was 
specified  as  being  separate  fi'om  the  wall  fire  in  the  input  file. 

real  GetFluxPI(ycGF)  -  Calculates  the  burner  contribution  to  the  surface  incident  heat  flux  if  it  was  specified 
that  the  burner  use  a  different  heat  flux  correlation  than  the  wall  fire  in  the  input  file. 

realyc(integer  nj)  -  calculates  the  absolute  y  coordinate  for  element  (x,nj). 

real  xcfinteger  ni)  -  calculates  the  absolute  x  coordinate  for  element  (ni,y). 

3.1.3  Input  and  Output  Files 

projecLdat  -  cont^  input  data  specifying  the  analysis  type,  material  properties,  dimt»nginng  and  desired  heat 
flux  and  flame  height  correlations  as  well  as  time  data. 
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heat_ga5.dat  -  Contains  heat  of  gasification,  depth  pairs  calculated  from  Cone  Calorimeter  data. 

tempout  -  Contains  a  list  of  nine  elements  (specified  by  vertical  number)  for  which  the  temperature,  timfi 
history  is  desired  (used  only  in  debug  mode). 

heatreaLdat  -  Contains  output  heat  release,  time  values. 

Jlame.dat  -  Contains  output  flame  height,  pyrolysis  front  height,  burnout  front  height,  flame  height  above 
floor  versus  time  values. 

debug.dat  -  Contains  mass  loss  rate  versus  time  data  (output  only  in  debug  mode). 

flux.dat  -  Contains  elemental  net  and  incident  heat  flux  data  versus  time  (output  only  in  debug  mode). 

t^^s.dat  -  Contains  temperature  versus  time  data  for  nine  elements  listed  in  tempout  (output  only  in  debug 


4.0  RESULTS  OF  THE  MODEL 

The  model  results  were  compared  with  full-scale  tests  to  evaluate  the  capabilities  of  the  computer 
program  in  predicting  vertical  spread.  Heat  release  rates,  surface  temperatures,  and  flame  hwghts  were 
compared  with  data  available  in  the  literature.  No  tests  were  found  that  described  in  adequate  all  of 
me  potential  characteristics  for  comparison.  In  addition,  the  test  conditions  were  not  adequately  d'Knibfd 
m  terms  of  the  ignition  temperatiue,  material  properties,  and  heat  flux  conditions.  Comparisons  were 
using  values  or  data  from  other  literature  sources.  Comparisons  were  also  made  to  assess  the  sensitivity  of 
me  inodel  to  several  parameters,  including  material  properties  (thermal  inertia  (kpc),  heat  of  gasification, 
igmtion  and  pyrolysis  temperature),  heat  flux  and  flame  height  correlations.  It  was  demonstrated  that  the 
predictions  can  be  altered  significantly  using  parameters  vrithin  the  bounds  of  reasonable  values. 
Consequently,  a  successful  prediction  relies  on  accurately  assessing  these  parameters  for  each  case. 

Validation  of  the  numerical  solution  for  upward  flame  spread  was  made  by  comparison  with  full-scale 
experimratal  results  using  material  properties  from  bench-scale  data  available  in  the  literature.  Full-scale  data 
was  available  for  v^ical  flame  spread  on  polymethylmethacrylate  (PMMA),  plywood,  and  wood  particle 
board.  The  present  numerical  solution  was  compared  to  the  experimental  measurement  of  flame  spread  over 
a  vertical  PMMA  surface  by  Wu,  Delichatsios,  and  de  Ris  (reference  (11)),  the  experimental  measurements 
of  flame  spread  over  a  vertical  surface  of  plywood  with  externally  applied  radiation  flux  by  Delichatsios  et 
al.  (reference  (12)),  and  vertical  flame  spread  measurements  on  wood  particle  board  by  Saito,  Quintiere  and 
Williams  (reference  (13)).  ’ 

4.1  Model  Inputs 

The  i^ition,  flame  spread  properties  and  dimensions  of  PMMA,  plywood,  and  wood  particle  board 
are  tabulated  in  Table  1.  The  properties  of  all  materials  in  Table  1  are  based  on  the  bench-scale  test  data 
available  in  the  literature.  The  sources  of  the  data  and  methods  used  for  this  determination  are  discussed  in 
the  subsequent  sections. 
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Table  1 .  Properties  and  Dimensions  of  Selected  Material 
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Data  reievaot  to  the  prediction  of  flame  spread  on  a  material  has  been  derived  experimentally  by 
several  investigators.  Experimental  results  by  Kashiwagi  and  OhlemiUer  (reference  (14))  for  the  radiant 
heating  of  PMMA  (Rohm  and  Hass  Type  G)  in  a  nitrogen  atmosphere  leads  to  a  leveling  of  the  surface 
temperature  as  vaporization  occurs.  For  PMMA,  these  temperatures  ranged  from  300-330°C  for  a  heat  flux 
of  17  kW/m^  and  365-380°C  for  40  kW/ml 

Tcwarson  and  Pion  (reference  (15))  and  Jackson  (reference  (16))  conducted  mass  loss  experiments 
for  estimating  the  steacfy  or  time-averaged  burning  rate  in  terms  of  heat  of  gasififiation.  Ah  for  PMMA 
(Rohm  and  Hass  Type  G).  The  heat  of  gasification  corresponding  to  the  Jackson  (reference  (16))  1600 
kJ/kg,  is  used  in  this  model. 

The  density  of  PMMA  has  been  measured  by  Quintiere  and  Harideroad  (reference  (10)),  and  Agrawal 
md  Atreya  (reference  (17))  as  well  as  by  Tewarson  and  Ogden  (reference  (18))  and  ranged  from  1 190-1200 
kg/m .  Tewarson  (reference  (19))  has  measured  flammability  parameters  of  materials  associated  with 
combustion  and  flame  spread.  In  the  experiments,  the  heat  of  combustion,  Ah,„  was  measured  to  be  24,900 
kJ/kg  for  PMMA.  These  properties  were  used  to  make  the  predictions  for  upward  flame  spread  for  PMMA. 

inertia  (kpc)  of  PMMA  was  found  by  fitting  ignition  data  of  Quintiere  and  Harkleroad 
Reference  (10))  using  an  ignition  temperature  of  320°C  (reference  (13)).  The  calculated  value  of  kpc  from 
^tiere  and  Harkleroad  (reference  (10))  ignition  data  was  0.60  (kW/m^.K)^  sec,  which  is  in  agreement  with 
the  value  of  0.57  (kW/m^-K)^  sec  reported  by  Tewarson  and  Ogden  (reference  (18)),  0.61  (kW/iti  -Kj 
^  obtained  by  Delichatsios,  Panagiotou,  and  Kiley  (reference  (20)),  and  0.61  (kW/m^-K)^  sec  reported  by 
Orloff,  de  Ris,  and  Markstein  (reference  (21)).  There  were  few  data  available  on  the  material’s  pyrolysis 

^perature.  A  value  50°C  -  60'’C  greater  than  the  ignition  temperature  is  representative  of  the  available 
data. 


The  ignition  and  flame  spread  properties  for  plywood  and  wood  particle  board  were  fojm 

the  experimental  bench-scale  data  of  Quintiere  and  Harkleroad  (reference  (10)),  Janssens  (reference  (22)), 
Parker  (reference  (23))  and  Quintiere  (reference  (7)).  The  thermal  inertia  was  calculated  by  fitting  critical 
Ignition  heat  flux  data  as  was  done  with  the  PMMA.  Fig.  9  and  Table  2  compare  the  predicted  and  actual 
igmtion  times  (from  Quintiere  and  Harkleroad,  1985  (reference  (10))  using  an  ignition  temperature  of  350°C 
as  reported  by  Janssens,  1991  (reference  (22))  and  a  thermal  inertia  of  0.475  (kW/m^-K)*  sec. 


Table  2.  Comparison  of  Experimental  and  Predicted  Ignition  Time  for  Plywood  (reference  (10)) 


kpc  =  0.475  (kW/m*-K)*  sec  T^  =  350  (®C) 

(Deduced  Value)  Tp  =  427  (®C) 

Janssens,  1991  (Reference  f22'0 

Erqrerimental  heat  uux  exposure 
(kW/m*) 

Measured  ignition  time 
(sec) 

Predicted  ignition  time 
(sec) 

17 

308 

310 

20 

185 

189 

24 

115 

114 

30 

62 

65 

41 

30 

31 

..  21 

20 

21 


Hgure  9  Incident  heat  flirc  versus  measured  and  predicted  ignition  time  for  plywood  at 
kpc  =  0.475  (kW/nf  K)^  sec  and  Tjg  =  350  °C. 


The  heat  of  gasification  for  plywood  was  deduced  from  the  experimental  results  of  Janssens 
(reference  (22))  and  Parker  (reference  (23)).  This  was  done  by  holding  all  material  properties  (  Table  1  and 
Appendix  ^  for  plywood)  constant  and  varying  the  heat  of  gasification  until  all  of  the  data  were 
satisfactorily.  The  pyrolysis  temperature  for  wood  was  given  in  Janssens  (reference  (22))  as  427“C.  Table 
3  compares  the  measured  and  predicted  heat  release  rate  with  Ahg  =  6850  kJ/kg  at  several  external  flux  levels. 


Table  3.  Comparison  of  experimental  and  predicted  heat  release  rate  results  for  plywood  at  Ah.  =  6850 
kJ/kg  (deduced  vdue) 


E^qrerimental  heat  flux 
e?qx>sure 
(kW/m") 

Average  heat  release  rate 
(kW/m") 

Predicted  heat  release  rate 
(kW/m^ 

Ahg  =  6850  kJ/kg 
(Deduced  Value) 

25  [22] 

85.20 

84 

35  [22] 

104.10 

104 

50  [22] 

142.10 

135 

50  [23] 

124.0 

135 

60  [22] 

156.30 

155 

Wood  particleboard  was  assumed  to  have  the  same  material  properties  as  plywood  with  the  exception 
of  Ae  thermal  inertia.  As  with  the  PMMA  and  the  plywood,  this  value  was  estimated  by  fitting  the  data  of 
Quintiere  (reference  (7)).  It  was  determined  that  the  same  thermal  inertia  (kpc  =  0.475  (kW/m^-K)*  sec) 
^timated  fen-  plywood  is  also  appropriate  for  particleboard.  A  comparison  of  the  measured  and  predicted 
i^ition  times  for  wood  particleboard  vrith  kpc  =  0.475  [kW/m^-K]^  sec  and  Tig  =  350®C  is  given  in  Table 
4. 
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Table  4.  Comparison  of  Experimental  and  Predicted  Ignition  Time  for  Wood  Particleboard 

(reference  (7)) 


kpc  =  0.475  (kW/m^-K)2  sec  T^=  350°C 

(Deduced  Value)  Tp  =  427°C 

(reference  (13)) 

Heat  Flux  (kW/m^) 

Measured  ignition  time  (sec) 

Predicted  ignition  time  (sec) 

20 

245-258 

399 

25 

155-180 

199 

30 

80-110 

120 

36 

68 

75 

40 

60 

58 

45 

49 

45 

50 

40 

35 

58 

30 

26 

60 

31 

24 

65 

18 

20 

4.2  Comparison  of  Predicted  with  Experimental  Results  for  Non-charring  Polymethylmethacrylate 

(PMMA) 

Wu,  Delichatsios,  and  de  Ris  (reference  (1 1))  conducted  experiments  to  measure  heat  release  rate  for 
bendi-scale  and  full-scale  PMMA  wall  fires.  In  the  bench-scale  tests,  a  vertical  PMMA  slab  measuring  0.90  m 
high,  0.2  m  wide  and  25  mm  thick  was  ignited  at  the  bottom  of  the  panel  by  20  ml  of  methanol  and  cotton  balls 
in  an  aluminum  dish  (0.025  m  x  0.20  m  x  0.01  m  high).  The  full-scale  flame  spread  experiments  were  carried 
out  with  a  25  mm  thick  PMMA  slab,  0.58  m  wide  x  5  m  high.  The  ignition  source  was  35  ml  of  heptane  in  a 
cqipo'  dish  (0.025  m  x  0.6  m  x  0.025  m  high)  at  the  bottom  of  the  wall.  This  was  simulated  in  the  model  as  a 
1  kW/m  line  fire  source. 
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Fig.  10  shows  a  comparison  of  heat  release  rate  during  upward  flame  spread  over  0.90  m  x  0.20  m 
vertical  PMMA  surface  (refer  to  Appendix  A,  Table  A-1  for  the  input  data).  The  solid  line  is  the  model 
predictiai  while  the  dashed  line  is  the  experimental  data  of  Wu,  Delichatsios,  and  de  Ris  (reference  (11)).  The 
properties  listed  in  Table  1  were  used  for  the  PMMA.  Due  to  the  sensitivity  of  the  computer  model  to  the  heat 
flux  and  flame  hei^  correlations,  sevaal  simulatitMis  of  the  test  were  perform  in  order  to  determine  the  optimum 
combination.  The  results  in  Fig.  6  show  reasonable  agreement  with  experimental  data  using  the  flame  height 
conelation  given  by  yj  =  0.052  Q'^  [Delichatsios  (reference  (4))],  and  the  flame  heat  flux  correlation  given  by 
Equations  (7  )  &  (8)  in  Section  2.2.  The  flame  heat  flux  in  the  burning  region  was  2 1  kW/m^.  Exact  values  of 
flame  heat  flux  after  igniticm  were  not  described  in  the  test  conducted  by  Wu,  Delichatsios,  and  de  Ris  (reference 
(11)).  A  value  of  2 1  kW/m^  heat  flux  after  ignition  was  found  to  yield  the  best  fit.  It  should  be  noted  that  the 
heat  flux  from  the  flame  to  the  wall  in  the  burning  region  is  not  actually  constant  but  rises  to  some  peak  value 
then  declines. 

Fig.  1 1  compares  the  predicted  heat  release  rate  and  the  measured  heat  release  rate  for  a  5.0  m  x  0.58 
m  votical  surface  of  PMMA  (Wu,  Delichatsios,  and  de  Ris  (reference  (1 1)).  The  input  file  is  listed  in  Appendix 
A.  The  solid  line  represoits  the  model  prediction  and  the  dashed  line  shows  the  experimental  results.  The  flame 
height  correlation  used  to  model  the  full  scale  test  was  yj  =  0.0666  Q'^,  developed  by  Tu,  and  Quintiere 
(referoice  (5)).  The  flame  height  correlation  gives  a  slightly  larger  flame  height  than  the  one  used  to  model  the 
bench-scale  tests  and  indicates  that  the  size  of  the  test  may  influence  the  height  of  the  flame  for  a  given  heat 
release  rate  per  unit  area.  The  heat  flux  correlation  used  in  the  full-scale  test  predictions  is  given  by  Equations 
(7)  and  (8).  A  3 1  kW/m^  heat  flux  was  used  after  the  material  ignited.  As  with  the  flame  height  correlation,  the 
heat  flux  was  greater  than  in  the  bench-scale  tests. 

4.3  Comparison  of  Predicted  and  Experimental  Results  for  Plywood 

Delidiatsios  et  al.  (reference  (12))  have  axiducted  a  theoretical  and  experimental  analysis  of  upward  fire 
spread  along  2.4  m  high,  0.61  m  wide  and  12.7  mm  thick  vertical  plywood  wall.  Five  full-scale  flame  spread 
e?q)erim^ts  were  ccaiducted.  The  samples  were  e?qx>sed  to  a  specified  heat  flux  from  a  large-scale  radiant  panel. 
The  following  measurements  were  made: 

(1)  Rate  of  heat  release, 

(2)  Total  heat  flux  to  the  specimen  surface, 

(3)  Surface  temperature,  and 

(4)  Propagation  of  the  pyrolysis  firont. 

The  fiill-scale  test  results  demonstrated  the  sensitivity  of  the  flame  spread  rate  to  the  external  heat  flux. 
The  samples  of  plywood  in  the  experiments  were  ignited  by  a  red-hot  nichrome  wire.  The  wire  was  preheated 
using  a  welder  power  supply  set  to  provide  a  40  Amp  current  through  the  wire.  After  preheating,  the  wire  was 
brou^t  into  direct  ccxitact  with  the  specimen  at  a  location  approximately  25  mm  above  the  base  of  the  specimen. 
A  thin  spacer  located  between  specimen  holder  and  the  center  of  the  specimen  was  used  to  make  the  specimen 
slightly  convex  at  the  base  and  ensure  good  contact  for  the  ignition  wire  over  the  entire  width  of  the  specimen. 
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Figure  10  Comparison  of  heat  release  rate  predictions  for  a  0.90  m  x  0.20  m  vertical 
PMMA  surface  with  Wu,  Delichatsios,  and  de  Ris,  1993 
(reference  (11)) 
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Time  (sec) 

Figure  1 1  Comparison  of  heat  release  rate  predictions  for  a  5.0  m  x  0.58  m 
PMMA  sur&ce  with  Wu,  Delichatsios,  and  de  Ris,  1993  Hats 
(reference  (11)) 


Comparisons  of  the  upward  flame  spread  model  to  the  full>scale  experimental  results  (Delichatsios 
ct  al.  (reference  (12))  are  shown  in  Fig.  12  and  14.  The  properties  of  pl)wood  used  for  the  prediction  were 
discuss^  in  Section  4.1.  The  heat  release  rate  during  pyrolysis  was  predicted  and  compared  to  the  measured 
results  in  Fig.  12,  (refer  Appendix  A,  Table  A-3,  test  1  for  input  data)  for  good-one-side  (GIS)  plywood.  In 
Fig.  12,  the  agreement  wiA  the  experimental  measurement  is  very  go^.  Typical  measmed  peak  heat  release 
rate  r^ts  in  this  test  series,  along  with  the  predicted  peak  heat  release  rate  results,  are  summarized  in  Table 
5.  Fig.  13  shows  the  heat  release  rate  curves  obtain^  from  the  computer  model  for  each  type  of  plywood. 


Table  5.  Comparison  of  Experimental  and  Peak  Heat  Release  Rate  Results  (reference  ( 1 2)) 


Test# 

Plywood 

Type 

Heat  Flux 
Exposure 
(kW/m*) 

Peak  Heal 
Flux 
(kW/m=) 

Initial 

Surface 

Temperature 

CO 

Measured  Peak 
Heat  Release 
Rate(kW) 

Predicted  Peak 
Heat  Release 
Rate(kW) 

1 

GIS 

5.2 

34-40 

160 

118 

114 

2 

RG 

11 

45-50 

250 

250 

136 

3 

RG 

mm 

40-45 

205 

130 

124 

4 

GIS 

4.8 

33-38 

155 

105 

100 

5 

GIS 

7 

40-45 

200 

175 

123 

GIS  -  Good-one-side 
RG  -  Rough  grade 


The  predicted  and  experimental  propagation  of  the  pyrolysis  front  (determined  using  surface 
thermocouples)  are  shown  in  Fig.  14.  In  general,  there  is  a  reasonable  agreement  between  predicted  and 
measured  propagation  of  the  pyrolysis  zone. 

Fig.  15  shows  the  temperature  rise  history  of  the  specimen  with  a  uniform  external  flux  of  5.2  kW/m^ 
over  a  penod  of  1200  sec  (20  min).  The  temperature  rise  and  steady  preheat  temperature  were  well  predicted 
by  the  model. 
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Heat  Release  Rate  (kW) 


Figure  13  Predicted  heat  release  rate  for  a  2.4  m  x  0.61  m  vertical  plywood  surface 
(reference  (12)) 


Height  (m) 


Time  (sec) 

Figure  1 4  Comparsion  of  pyrolysis  front  propj^ation  predictions  for  a  2.4  m  x  0  6 1  m 
vertical  plywood  surface  with  Delichatsios  et  al.  1994  data  (refwence  (12)) 
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Fig.  15  Comparison  ofprebeating  temperature  prediction  for  a  2.4  mx  0.61m  vertical 

pfywood  sur£me  with  Delichatsios  et  aL  1994  data  (reference  (12)) 


The  total  flame  heat  flux  after  ignition  was  measured  by  Delichatsios  et  al.(reference  (12))  for  each 
experiment.  For  the  prediction  of  heat  release  rate  for  each  specific  test  in  Figs.  12  -  14,  the  heat  flux  after 
ignition  is  calculated  as  average  peak  heat  flux  minus  imposed  heat  flux  from  Table  5.  Peak  heat  fluxes  after 
igniticm  fcx  all  tests  conducted  Ity  Delichatsios  et  al.  (reference  (12))  were  greater  than  heat  flux  data  summarized 
by  Quintiere  (reference  (7))  (refer  Section  2.2  Fig.  3).  Based  on  Fig.  3,  five  more  heat  flux  correlations  have 
been  developed  and  are  shown  in  Fig.  16,  with  correlations  1, 2,  and  3  discussed  in  Section  2.2.  Correlation  6, 
7,  8,  9,  and  10  in  Fig.  16  all  lie  above  the  heat  flux  plateau  in  Fig.  3.  Table  6  shows  heat  flux  after  ignition 
calculated  from  Table  5  with  heat  flux  correlation  used  by  the  model  to  compare  heat  release  rate  for  each  test 
coiducted  by  Delichatsios  et  al.  (reference  (12)).  The  flame  height  correlation  given  by  Equation  (2)  in  Section 
2.2  with  Cf  =  0.052  was  used. 


Table  6.  Heat  Flux  after  Ignition  and  Heat  Flux  Correlation  for  Plywood 


Test# 

Imposed  Heat  Flux 
(kW/m^) 

Heat  Flux  after  Ignition 
(kW/m^) 

Heat  Flux  Correlation 

1 

5.2 

37 

8 

2 

11 

38 

8 

3 

7.5 

38 

8 

4 

4.8 

33 

9 

5 

7.0 

38 

8 

4.4  Comparison  of  Predicted  and  Experimental  Results  for  Particle  Board 

Saito,  Quintiae,  and  Williams  (reference  (13))  have  examined  the  mechanisms  and  rate  of  upward  flame 
spread  along  thermally  thick  vertical  sheets  for  both  charring  and  non-charring  fuels.  Measurements  of  spread 
rates,  flame  heights  and  of  surface  temperature  histories  were  reported  for  PMMA  and  Douglas-fir  particle  board 
for  flames  initiated  and  siqjpcxted  by  a  line  source  gas  burner  located  at  the  base  of  the  panels.  The  heat  release 
rate  of  the  burner  was  varied  fixxn  8  kW/m  to  about  100  kW/m.  The  pyrolysis  front  was  determined  from  surface 
temperature  measuianaits  along  the  sanqjle  vertical  centerline  and  flame  heights  were  determined  visually.  The 
sustained  flame  spread  was  not  observed  for  wood  particle  board  samples  even  if  the  burner  was  left  on  for  ten 
minutes  in  the  experiments  conducted  by  Saito,  Quintiere,  and  Williams  (reference  (13)).  The  maximum  height 
of  the  char  yield  (pyrolysis  front)  increased  appreciably  with  increasing  burner  energy  supply. 

The  properties  of  wood  particle  board  (1.8  m  high,  0.30  m  wide  and  13  mm  thick)  discussed  in  Section 
4. 1  (refer  to  Table  1  and  Appendix  A)  were  used  to  make  the  predictions  shown  in  Fig.  17  and  Table  7  for 
various  intensities  of  the  line  burner. 


33 


45 


y/Xf 


Figure  16  Heat  flux  correlation  for  the  model 
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%  Menncd  pyrolysis  height 

Q  Predicted  pyrolysis  height 

V,  -  0.35,  kpc  -  0.475  (kW/m’  K)’  sec,  T^-  350  *0 

^  Predicted  pyrolysis  height 

V,-0.85,kpc-0.475  (kW/in’K)’sec,T^-350*C 


1.0  k 


Predicted  pyrolysis  height 
y,-0.0«I.Q“ 

T„-350*C 

T,-427*C 

kpc  -  0.475  (kW  /  K)*  sec 

q'-30kW/ni* 
dh  -6850ki/kg 


Figure  17  Comparison  of  pyrolysis  height  predictions  for  a  1.8  m  x  0.3  m  vertical 
wood  particle  board  with  Saito,  Quintiere,  and  Williams,  1985  data . 
(reference  (13)) 
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Ignition  by  the  line  burner  at  the  base  of  the  particle  board  sample  at  various  energy  releases  was 
simulated  by  a  ui^orm  constant  distribution  of  heat  flux  to  the  surface  over  the  flame  height.  The  flamw 
height  correlation  used  to  compare  the  prediction  with  data  is  directly  obtained  from  the  reported  by 
Smto,  Quintiere,  and  Williams  (reference  (13)).  Their  results  show  a  value  of  Cf  =  0.0421  m^/kW^  agree 
with  a  2/3  power  dependence  of  the  flame  height  on  the  heat  release  rate  per  unit  length.  The  flame  hftipht 
is  as  follows: 


2^ 

3y=  0.0421 


(24) 


The  predicted  maximum  pyrolysis  heights  using  a  char  yield  (7)  of  0.35  (reference  (24))  do  not  agree 
with  the  experimental  data  quantitatively.  Thqr  do,  however,  have  a  similar  trend  (slope  in  the  ^ta). 
the  char  yield  to  0.85  produces  results  in  better  agreement  with  the  experimental  data  of  Saito,  Quintiere,  and 
Williams  (reference  (13)).  This  indicates  that,  although  there  is  some  consistent  reproducible  aspect  to  the 
experimental  data,  one  or  more  phenomena  are  not  actually  being  predicted  and  are  being  compensated 
through  an  unrealistically  high  char  yield. 


Table  7.  Comparison  of  Experimental  and  Predicted  Pyrolysis  Height  Results  for  Wood  Particle  Board 
(reference  13)) 


Rate  of  heat 
release 
(line  fire  gas 
burner) 
(kW/m) 

Measured 
maximiun 
pyrolysis  height 
(7) 

(m) 

Predicted  maximum 
pyrolysis  height  (7) 

(m) 

7=0.35,Ti,=  350‘’C 
kpc  =  0.475  (kW/m^lQ^sec 

Predicted  maximum 
pyrolysis  height  (7) 

(m) 

7=0.85,T^  =  350°C 
kpc  =  0.475  (kW/m2.K)2sec 

22 

0.40 

0.65 

0.42 

38 

0.60  -  0.61 

0.95 

0.61 

58 

0.91-0.98 

1.25 

0.82 

76 

0.91 

1.48 

1.0 

96 

1.3 

1.75 

1.2 

4,5  Sensitivity  Analysis 

The  rate  of  flame  spread  over  a  material  surface  is  influenced  by  a  number  of  parameters  as 
in  preceding  sections.  Some  of  these  parameters  have  a  greater  impact  than  others,  and  consequently  more 
attention  should  be  given  to  these  when  modeling  a  particular  situation.  This  section  systematically  explores 
the  impact  on  the  model  results  for  seven  of  the  input  parameters; 


(1)  Heat  of  gasification, 

(2)  Thermal  inertia. 
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(3)  Ignition  temperature, 

(4)  Pyrolysis  temperature, 

(5)  Flame  height  correlation, 

(6)  Flame  heat  flux  correlation,  and 

(7)  Flame  heat  flux  after  ignition. 

pje  baseline  case  selected  was  the  small-scale  PMMA  predicted  results  (Fig.  10).  All  parameters 

^  being  examined.  The  sensitivity  to  the  seven  parameters  is 


The  eff^  of  variations  m  the  heat  of  gasification  are  shown  in  Fig.  18  in  which  predictions  for  heat 
of  gasification  ^m  the  actual  value  of  1600  kJ/kg  to  2000  kJ/kg  are  shown.  As  the  Figure  illustrates 
vanations  of  ^  relatively  smaU  magnitude  have  a  significant  effect  on  the  results.  This  points  to  the  need 
lor  excellent  characterization  of  the  material  properties  to  allow  accurate  predictions. 


n,xir,  2^1  vanations  m  the  thermal  inertia  from  the  actual  value  of  0.6  (kW/m^K)*  to  0  84 

(kW/m  K)  are  shown  m  Fig.  19.  The  sensitivity  of  the  results  to  variations  in  the  thermal  inertia  is  shown 
to  be  less  ^  that  associated  with  the  heat  of  gasification.  As  the  discussion  of  the  inputs  in  Section  4.1 
mdicated,  the  heating  and  ignition  behavior  of  materials  can  typically  be  well  represented  by  the  integral  heat 
transfer  model  and  the  igniUon  temperature  concept.  As  such,  significant  uncertainties  in  the  ignition 
respoMe  as  represented  by  the  thermal  inertia  and  ignition  temperature  are  not  expected.  Figure  20  shows 
®  var^g  the  ignition  temperature  input.  It  is  possible  to  determine  the  ignition  temperature  to 
witnm  20  C  so  that  significant  errors  due  to  ignition  temperature  uncertainties  are  not  large.  It  needs  to  be 
remernbered  that  the  thermal  mertia  and  ignition  temperature  are  used  together  to  predict  the  ignition  of  each 
numencal  elcmrat.  Thus,  the  real  key  to  assessing  the  uncertainties  is  the  degree  of  fit  in  the  time  to  ignition 
over  a  rruige  of  heat  fluxes  as  shown  in  Fig.  9.  ^ 

The  effects  of  pyrolysis  temperature  on  the  heat  release  rate  are  shown  in  Fig.  21.  The  pyrolysis 
pni^ly  governs  the  radiative  heat  losses  from  the  burning  surface.  For  some  materials  like 
PMMA,  the  surf^  temperature  during  burning  is  quite  constant.  Other  materials  like  wood  can  experience 
vi^ations  in  s^ace  temperature  dunng  burning  in  excess  of  100°C.  As  such,  the  uncertainties  associated 
wnth  toe  use  of  a  constant  pyrolysis  temperature  can  vary  widely  depending  on  the  material.  This  variable 
may  be  responsible  for  the  difficulties  in  predicting  the  flame  spread  distances  on  toe  wood  particleboard 


The  effect  of  the  flame  height  correlation  constant  is  demonstrated  in  Fig.  22.  The  range  of  constants 
exammed  m  Fig.  22  corresponds  to  the  range  of  constants  reported  in  the  literature.  As  toe  Figure  shows  this 
instant  can  ^edly  affect  the  results  of  the  prediction.  Additional  detailed  investigation  of  wall  flame 
neight  correlations  is  mdicated  by  this  sensitivity. 
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Hgure  1 8  Comparison  of  heat  release  rate  predictions  for  a  0.90  x  0.20  m  vertical 
PMIi^  sur&ce  at  different  values  of  heat  of  gasification  with  Wu, 
Delichatsios,  and  de  Ris,  1993  data  (reference  (1 1)) 
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Figure  1 9  Comparison  of  heat  release  rate  predictions  for  a  0.90  m  x  0.2  m  vertical 
PMMA  surface  at  different  values  of  thermal  inertia  with  Wu, 
Delichatsios,  and  de  Ris,  1993  data  (reference  (11)) 


Figure  20  Comparison  of  heat  release  rate  predictions  for  a  0.90  m  x  0.20  m  vertical 
PMMA  surface  at  different  ignition  temperature  with  Wu,  Delichatsios, 
and  de  Ris,  1 993  data  (reference  (1 1)) 
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Figure  21  Co^arison  of  heat  release  rate  predictions  for  a  0.90  m  x  0.20  m  vertical 
PM^  surf^  at  different  pyrolysis  temperature  with  Wu,  Delichatsios, 
and  de  Ris,  1993  data  (reference  (11)) 
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Figure  22  Comparison  of  heat  release  rale  predictions  for  a  0.90  m  x  0.20  m  vertical 
PMMA  surface  different  flame  height  correlation  with  Wu,  Delichatsios, 
and  de  Ris,  1993  data  (reference  (1 1)) 
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Figure  23  Comparison  of  release  rate  predictions  for  a  0.90  m  x  0.20  m  vertical 
PMMA  surface  different  heat  flux  correlations  with  Wu,  Delichatsios, 
and  de  Ris,  1993  data  (reference  (1 1)) 
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Figure  24  C^P^on  of  release  rate  predictions  for  a  0.90  m  x  0.20  m  vertical 


The  effect  of  variations  in  the  heat  flux  to  the  wall  during  burning  are  illustrated  in  Fig.  23.  The 
TMge  of  heat  fluxes  examined  in  the  Figure  corresponds  to  the  range  of  experimental  data  available  in  the 
literature.  As  the  Figure  indicates,  the  variations  in  heat  release  rate  are  very  significant.  This  also 
that  additional  work  in  the  measurement  and  correlation  of  heat  fluxes  to  walls  is  needed.  The  level  of 
certainty  represented  by  the  literature  is  simply  not  suflicient  to  support  accurate  predictions. 

The  effwt  of  heat  flux  to  the  alreacfy  burning  surface  is  shown  in  Figure  24.  The  range  of  fluxes 
shown  in  the  Figure  represents  the  range  of  values  reported  in  the  literature.  Clearly,  the  behavior  of  the 
material  is  critically  impacted  by  the  heat  flux  to  the  burning  material.  The  results  are  particularly  sensitive 
to  the  heat  flux  for  lower  heat  fluxes  where  spread  rates  are  slow  and  indefinite  flame  spread  is  not  erqiected. 
It  is  known  that  the  heat  flux  to  the  burning  surface  is  scale  dependent  with  larger  fluxes  vrith  taller 

wall  flames.  Howwer,  the  details  of  the  scale  effects  are  not  sufficiently  well  known  to  allow  their  inclusion 
in  the  model.  Additional  work  in  this  area  is  needed. 

4.6  Discussion  of  the  Sensitivity  Analysis 

The  inputs  ^sessed  in  the  prior  section  can  be  classified  generally  as  material  properties  inputs  and 
wall  heat  transfer  inputs.  It  can  be  hypothesized  that  the  wall  heat  transfer  inputs  are  only  modestly 
dependent  on  the  fire  properties  of  the  material,  though  the  luminosity  of  the  flame  is  clearly  dependent  on 
the  nature  of  the  fuel  volatiles.  The  range  of  wall  heat  transfer  inputs  examined  in  the  prior  section  generally 
reflects  the  range  of  results  available  in  the  research  literature  at  this  time.  While  the  use  of  average  values 
for  thwe  inputs  were  gmerally  quite  successfully  used  for  data  comparisons  in  Sections  4.2, 4.3,  and  4.4,  the 
sensitivity  analysis  indicates  that  improved  precision  in  these  inputs  is  needed  to  support  accurate  predictions 
of  flame  spread  and  resulting  heat  release  rate.  The  need  for  accurate  inputs  is  particularly  keen  under 
conditions  were  acceleratory  flame  spread  is  not  present  or  where  it  is  marginally  present.  Since  a  majOT  goal 
of  these  models  is  to  assess  materials  to  allow  selection  of  materials  which  do  not  support  acceleratory  flame 
spread,  this  need  for  increased  precision  is  critically  important  for  practical  use  of  the  and  is  in  no 
sense  academic.  This  points  to  the  need  for  well  controlled  experimental  woric  in  this  area. 

The  material  properties  inputs  are  equally  critical  to  the  performance  of  the  model.  As  Section  4.4 

indicates,  the  pcxirest  prechetions  of  material  performance  were  for  experimental  ennHitinn^  i^riiere 
acceleratory  flame  spread  did  not  occur  for  the  wood  particleboard.  While  it  was  anriripateH  that  the 
relatively  simple  solid  phase  material  response  models  used  would  be  sufficiently  accurate  to  allow  good 
predictions  of  material  behavior,  it  is  clear  at  this  point  that  more  detailed  models  will  be  required  to 
reproduce  actual  inaterial  fire  performance.  The  ignition  model  seems  to  be  fully  adequate  at  this  time. 
However,  the  burning  rate  model  used  here  needs  to  be  further  developed  to  achieve  the  required  accuracy. 
The  pyrolysis  tempCTature  assumed  to  be  constant  in  this  model  can  vary  enough  to  seriously  compromise 
the  restdts  for  materials  like  wood.  Further,  the  heat  of  gasification  model  also  needs  improvement  to  achieve 
the  desired  level  of  accuracy. 

The  difficulties  with  using  a  more  sophisticated  model  of  the  solid  phase  response  of  the  burning 
material  are  two-fold.  First,  as  the  models  are  made  more  sophisticated,  they  also  become  less  general.  This 
means  that  there  may  need  to  be  multiple  models  applicable  to  different  material  types.  For  instance  a 
septate  rnodel  may  be  required  for  charring  and  non-charring  materials  or  thick  and  thin  fuels.  This  is  an 
obvious  difficulty  in  that  the  ease  of  use  of  the  model  can  quickly  be  degraded.  The  difficulty  is  the 
increas^  data  requirements  of  more  sophisticated  models.  Ultimately,  this  concern  is  greater  than  the  model 
generality  con^.  As  more  inputs  are  required,  the  number  and  detail  of  the  bench  scale  tests  needed 
increases.  This  would  ultimately  require  many  more  tests  and  require  the  cone  calorimeter  mgpiimmtntion 
to  be  developed  to  a  much  higher  level  than  is  currently  used.  In  terms  of  the  ultimate  use  of  the  model,  the 
demands  on  the  number  and  sophistication  of  the  bench  scale  tests  required  are  the  most  serious  conceriL 
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While  the  foregoing  discussion  focuses  on  the  uncertainties  and  limitations  of  the  present  model,  the 
comparisons  shoAvn  in  Sections  4.2, 4.3,  and  4.4  arc  very  good  indeed.  The  use  of  the  best  available  inputs 
and  parameters  h^  been  shown  to  result  in  very  credible  predictions.  Thus,  we  have  direct  evidence  of  the 
value  of  the  modeling  approach  adopted  here  and  the  prospects  for  a  valuable  material  assKssm^t  tool. 


5.0  PREDICTION  OF  VERTICAL  FLAME  PROPAGATION  OVER  THERMALLY  THIN 
MATERIAL  ON  NON-COMBUSTIBLE  SUBSTRATES 

5.1  Introduction 


The  goal  of  the  flame  spread  model  is  to  provide  a  link  between  bench  scale  test  results  and  fiill  scale 
fire  performance.  The  current  model  has  only  been  developed  to  date  to  include  upward  flame  spread  on 
vertical  surfaces.  Preliminary  work  to  generalize  this  model  is  described  in  Section  5.4.  However,  the 
development  of  the  model  and  the  room  comer  fire  testing  performed  onboard  the  ex-USS  SHADWELL 
(reference  (26))  using  the  U.S  Navy  material  test  program  has  motivated  a  simple  method  for  correlating 
bench-scale  fire  test  results  with  full  scale  comer  fire  test  results.  The  correlational  method  described  here 
was  developed  previously  and  reported  in  NRL  Letter  Report  (reference  (27))  and  is  included  here  for 
completeness. 

The  flame  spread  performance  of  U.S.  Navy  materials,.  Williams,  et  al.  (reference  (28))  ,  (Glass 
Reinforced  Plastic  Nomex  Panel,  Manville  Thermal  Insulation,  Imi-Tech  Acoustic  Insulation,  and  Waffle- 
Board  Acoustic  Insulation)  have  been  evaluated  and  compared  with  the  flame  spread  test  results  on  textile 
wall  covering  materials  on  non-combustible  substrates  (Harkleroad,  reference  (29)).  The  Navy  cone 
calorimeter  tests  (reference  (28))  were  conducted  by  the  Carderock  Division,  Naval  Surface  Warfare  Center 
(NSWC/CD)  Code  643,  and  full-scale  calorimeter  tests  were  conducted  aboard  ex-USS  SHADWELL  (LSD- 


In  the  next  section,  the  flame  spread  model  is  presented  to  provide  the  theoretical  basis  for  the 
flammability  parameter  developed  by  Mowrer  and  Williamson  (reference  (25)). 

5J2  Formation  of  Flame  Spread  Parameter 

The  flame  spread  model  developed  by  Mowrer  and  Williamson  (reference  (25))  is  based  on  the 
approach  presented  by  Quintiere,  Harkleroad  and  Hasemi  (reference  (30))  for  the  prediction  of  upward  flame 
spread  ovct  materials.  The  upward  flame  spread  model  (reference  (25))  includes  consideration  of  the  finite 
burning  time  tj,  of  thin  fitels.  The  consumption  of  the  wall  fuel  results  in  burnout  of  the  flame  at  that 
location,  wWch  is  an  important  aspect  of  flame  spread  on  thin  fuels. 
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In  this  model,  the  flame  spread  rate  is  deflned  as  the  rate  of  pyrolysis  front  advance; 

V  =  ^  9  ~  =  yp^  -  yp) 

’’  dt  t,  ~  t. 


v^ere  Vp  is  the  upward  spread  velocity  (m/s),  ^  is  the  height  of  the  pyrolysis  front  (m),  t  is  time  (s)  t,  is  the 
fl^e  spread  tune  (s),  and  is  the  flame  height  (m).  The  characteristic  flame  spread  (or  ignition)  time  is 
aetmed  m  terms  of  a  simple  thermal  model  of  heating  a  wall  with  constant  thermal  properties 


/tpc 


il,  -  T) 


where  aU  terms  have  been  previously  defined.  Once  burnout  begins,  the  velocity  of  the  burnout  front  can  be 
expressed  as: 


V  -  yp  *  y,(t)  -  yp) 

‘‘  dt  t. 


wher^  Vb  is  the  upward  velocity  of  the  burnout  front,  t  bo  is  the  burning  duration(s),  t,  is  the  burnout 
time(s),  and  yb  is  the  burnout  front  height  (m). 

A  linea^flame  height  approximation  is  used  to  describe  the  flame  height  required  in  Equation  (1),  following 
iaito,  Quml^,  and  WUliams  (reference  (13)),  Quintiere,  Harkleroad,  and  Hasemi  (reference  (30)),  and 

eary  Qimtiere  (reference  (31)).  Before  burnout  begins,  this  flame  height  approximation  can  be 
expressed  as  follows: 

y,  V  ^  (yv\ 


After  burnout  begins,  the  dimensionless  flame  height  is  expressed  as  follows: 

(29) 

The  parameter  Iqis  a  correlating  factor  used  to  define  the  flame  length.  Cleary  and  Quintiere  (reference  (3 1)) 

mggest  a  value  of  approximately  0.01  m*/kW,  for  Iq.  Using  Equation  (28)  for  times  t  <  tb.  Equation  (25)  can 
be  rewritten  as  follows:  v  /vau 


ik^E" 


47 


(31) 


Equation  (30)  can  be  integrated,  with  the  limit  that  y  =  ypo  at  t  =  0  and  y  “  yp  at  t: 


yp  =  ypo  exp 


-  1)^ 


Equations  (30)  and  (3 1),  together  with  Equation  (28),  suggest  that,  before  burnout,  the  flame  spread  rate 
will  accelerate  if  y,  >  yp  and  decelerate  if  yf  <  yp,  i.e.,  if  k,  E"  <  1 .  After  burnout  begins,  at  times  t  >  tj,,  the 
net  rate  of  flame  propagation  can  be  expressed  as  the  difference  of  the  pyrolysis  front  velocity  and  the 
burnout  front  velocity: 


K(0  = 


y^ 


9 


(32) 


Using  Equation  (29),  Equation  (32)  can  be  rearranged  to  the  fbllowing: 


^(yp  -yt)  =  O';,  -y^) 


{k,E"  -  1)  -  r; 


ho 


(33) 


Eqi^on  (33)  can  then  be  integrated,  with  the  limit  of  (yp  -  yj  =  (yp,  -  ypj  at  t  =  tt  and  typ  -  yO  =  typ  - 
at  time  t,  to  yield  the  pyrolysis  zone  height: 


Equation  (34)  suggests  that,  following  the  onset  of  fuel  burnout,  the  potential  for  apr^lftrating  flame 
spread  depends  on  a  balance  among  three  parameters:  the  normalized  flame  height,  (yr  -  yb)  /  (y  *  y  b), 

^ch  is  represented  per  Equation  (29)  as  a  linear  function  of  unit  heat  release  rate ,  E;  the  flame  spr^ 
time,  tft  given  by  Equation  (26);  and  the  burning  duration,  tbo-  If  the  parameter  (kf  £"-  1,  the 

predicted  flame  spread  will  accelerate.  The  flammability  parameter  is  thus  defined  as  E"kf-t/tb. 

While  this  model  is  based  on  several  idealizations,  it  is  expected  that  this  flammability  parameter 
characterizes  a  material  with  regard  to  vertical  flame  spread.  However,  attention  mygt  be  paid  to  the 
methods  used  to  evaluate  E",  t^  and  tb.  Mowicr  and  Williamson  evaluated  E"  as  the  peak  heat  release  rate 
of  the  material,  tf  as  the  ignition  time,  and  tb  as  the  time  from  ignition  to  the  peak  heat  release  rate.  Thty 
evaluated  these  quantities  at  both  cone  calorimeter  heat  fluxes  (30  and  50  kW/m^  tested  by  Harkleroad 
(reference  (29)).  Based  on  their  results,  the  correlation  was  better  using  the  cone  data  at  50  kW/m^ 

There  are  problems  with  the  methods  proposed  by  Mowrer  and  Williamson  for  deducing  E"  and  tb 
from  cone  calorimeter  data.  The  problems  are  umque  to  thin  materials  which  have  bum  times  less  than 
approximately  a  minute.  The  time  response  characteristics  of  the  cone  calorimeter  are  such  that  the  peak 
measured  heat  release  rate  is  less  than  the  actual  peak  for  these  thin  materials  due  to  the  small  bum  timf!  The 
effect  of  various  burning  durations  can  be  seen  from  Fig.  25  for  a  methane  burner  at  a  heat  release  rate  of  6.80 
kW  operated  for  various  durations.  Of  course,  if  the  cone  had  a  zero  response  time,  the  measured  heat  release 
rates  would  be  square  wave  pulses  with  the  width  equal  to  the  burning  duration.  For  the  longer  burning 
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dwatio^  (tert  1, 120  sec  bum  duration)  the  actual  heat  release  rate  is  measured  after  about  10  seconds  For 
shorter  b^  du^ons  (test  5, 10  sec  and  test  6, 5  sec),  the  peak  recorded  heat  release  rate  occurs  at  about  5- 
10  s^nds,  and  the  actual  burning  rate  is  never  recorded.  While  the  response  time  of  the  gas  analysis  system 
on  the  cone  calorimeter  does  not  allow  correct  measurement  of  the  heat  release  rate,  there  is  hope  that  the 
cumulative  heat  release  may  be  measured  correctly  despite  the  time  response  limitations  of  the  system.  The 
cumulative  hwt  release  is  the  area  under  the  heat  release  rate  verses  time  curve.  Table  8  shows  the  predicted 
and  measured  cumulative  heat  releases  for  the  various  bum  durations.  The  results  indicate  that  the  cone 
calorimeter  can  correctly  measure  the  cumulative  heat  release  for  short  duration  bums. 


Table  8.  Measured  and  Predicted  Cumulative  Rate  of  Heat  Release  for  a  Methane  Burner 

Operated  at  6.8  kW 


Test  number 

Test  duration 
(sec) 

Measured  cumulative  heat 
release  rate  ■ 

(kJ) 

Precheted  cumulative  heat 
release  rate 

(kJ)  1 

1 

120 

830.73 

816 

2 

60 

429.81 

408 

3 

30 

204.89 

204 

4 

15 

109.47 

102 

5 

10 

77.75 

68 

6 

5 

36.70 
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problems  with  the  bum  time  used  by  Mowrer  and  Williamson,  which  is  taken  as  the 
time  from  Ignition  to  the  time  of  the  peak  rate  of  heat  release  rate.  There  are  several  problems  with  this 
memod.  The  first  IS  the  significant  effect  of  response  time  on  the  measured  time  to  peak  heat  release.  This 
IS  directly  related  to  the  peak  heat  release  measurement  concern  and  is  particularly  acute  since  peak  heat 
release  ratw  typicaUy  occur  shortly  after  ignition.  Second,  if  the  peak  heat  release  occurs  soon  after  ignition 
u  release  occurs  after  the  peak,  the  actual  heat  output  duration  may  not  be  well  represented 

by  the  Mowtct  ^  Williamson  definition  of  bum  time.  Typically,  thick  and  thin  coverings  of  the  same 
matenal  would  have  the  same  bum  time  as  determined  by  the  Mowrer  and  Williamson  method  whereas  their 
Observed  burning  durations  would  be  very  different. 

..  .  ^  to  address  these  problems,  new  methods  of  cone  calorimeter  data  reduction  have  been 
stacked,  and  a  mocMed  method  has  been  developed.  The  bum  time,  t^,  is  taken  as  the  time  fi-om  ignition  until 
the  matCTial  stops  flammg.  This  is  best  determined  visually,  but  can  be  determined  from  the  heat  release  rate 
verses  time  output  from  the  cone  calorimeter.  For  all  materials  examined  by  Harkleroad  (reference  (29)), 
burn  times  were  determined  from  the  cone  calorimeter  heat  release  rate  data.  For  the  Navy  materials 
tae  bum  time  was  determined  from  visual  observation  during  the  cone  calorimeter  test  The  heat  release  rate, 
E  ,  IS  taken  as  the  cumulative  heat  release  as  routinely  determined  in  the  cone  calorimeter  divided  by  the  burn 
tune,  tb.  This  is  an  average  heat  release  rate  for  the  material  during  the  active  burning  period  These  morfifipH 
methods  avoid  some  of  the  experimental  difficulties  with  the  cone  calorimeter  as  appUed  to  thin  materials. 
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5  J  Experimental  Data  and  Implication  of  the  Model 

The  flammability  parameter  derived  from  cone  calorimeter  test  results  for  the  50  kW/m^  heat  flux 
exposure  for  nmc  textile  wall  covering  materials  on  gypsum  board  are  tabulated  in  Table  9  along  with  the 
pe^  heat  rele^  rato  measured  m  the  large-scale  comer  tests  performed  by  Fisher,  MacCraken,  and 
Wilhamson  (referenw  (32))  at  the  Umversity  of  California,  Berkeley.  The  flammability  parameter  from  cone 

S ^  protection  (PFP)  Navy  materials  are  tabulated  in  Table  10  for 

50  kW/m  heat  flux  exposure  to  the  specimen  surface  along  with  the  peak  heat  release  rates  measured  in  the 
full-scale  rorncrtKts  performed  on  the  ex-USS  SHADWELL.  The  correlation  of  the  peak  full-scale  heat 
release  with  the  flammability  parameter  from  data  for  the  50  kW/m^  exposure  level  for  each  of  the  Navy 
matoals  and  the  mne  ^e  wall  covering  materials  on  gypsum  board  arc  shown  in  Fig.  26.  Fig.  26  shows 
the  correlation  of  the  full-scale  results  by  the  flammability  parameter  is  quite  good.  Based  on  the 
perfom^cc  required  m  full-scale  appUcations,  Fig.  26  could  be  used  to  establish  the  required  value  of  the 
tiammabihty  parameter  for  mclusion  m  material  specifications  for  that  application. 


Table  9.  Summary  of  Flammability  Parameter  for  Textile  Wall  Covering  Materials  on  Gypsum  Board 
50  kW/m  Incident  Heat  Flux  Exposure  (reference  (29),  (3 1)) 


*  -  See  Fig.  26  for  description 

**  -  Parenthetical  peak  heat  release  rates  are  for  two  foot  wide  samples  while  normal  values  are  for  one  foot  wide 
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Fig.  26  Coiiq)arison  of  Navy  results  with  textile  wall  covering  on  gypsum  board  results 
(Harideroad,  1989,  reference  (29)  and  Fisher  ct  al.  1986  (reference  (32)) 
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Table  10.  Summaiy  of  Flammability  Parameter  for  the  PFP  Navy  Test  Materials, 
50  kW/m^  Incident  Heat  Flux  Exposure  (reference  (28)) 


PFP  Navy 
test  material 

Heat 

release  rate 
from  full- 
scale  fire 
test 
(kW) 

Bsm 

Flame 

length 

parameter 

kf 

(m^/kW) 

Flame 
spread  time 
tf 

(s) 

Burnout 

time 

tb 

(S) 

Flammability 
Parameter  = 

[E"k,-tA] 

GRPNomex 

panel 

70 

68 

0.01 

14 

120 

0.56 

Manville 

thermal 

insulation 

80 

5 

0.01 

8 

20 

-0.35 

Imi-tech 

acoustic 

insulation 

20 

■ 

0.01 

9 

10 

-0.82 

WaflQe  board 
acoustic 
insulation 

150 

60 

0.01 

9 

20 

0.15 

If  one  examines  Fig.  26,  the  GRP  data  point  appears  to  be  somewhat  of  an  outlier.  Based  on  the 
flammability^  parameter  of  0.56,  one  would  have  expected  a  higher  full-scale  heat  release  rate.  This  material 
burns  vigorously  where  the  surface  is  breached  by  cracking  due  to  heating.  This  behavior  is  emphasized  in 
the  cone  calorimeter  in  that  such  a  crack  alw^s  occurred  in  the  small  sample  while  cracking  occurred  at  a 
relatively  ^der  spacial  separation  in  the  hiU-scale  test.  Thus,  the  cone  calorimeter  test  results  are  worse  than 
expected  in  full-scale.  Such  scaling  problems  must  be  expected  to  occur  for  some  materials  in  a  small-scale 
test  and  is  an  uAermt  limitation  of  small-scale  testing.  However,  the  effect  is  not  significant  for  any  of  the 
materials  tested  in  this  program 

A  flammability  parameter  based  on  concurrent  vertical  flame  spread  modeling  has  been  developed 
to  allow  correlation  of  bench-scale  ignition  and  heat  release  rate  measurements  in  the  cone  raiorim^tw  with 
full-scale  fire  p^ormance.  The  cone  calorimeter  data,  at  an  exposure  level  of  50  kW/m^  for  both  Navy 
insulation  materials  and  the  textile  wall  coverings  on  gypsum  board  shown  in  Figure  26,  demonstrate  the 
general  trend  that  materials  with  low  heat  release  in  full-scale  fire  tests  have  low  flammability  parameter 
values,  and  materials  with  high  heat  release  rate  in  real-scale  fire  tests  have  high  values  of  the  flammability 
parameter.  As  expected,  the  behavior  of  materials  change  at  a  flammability  parameter  of  about  one.  The 
®®^^l2tion  successfully  luiks  cone  calorimeter  results  to  full-scale  fire  performance  for  Navy  wall  covering 
material,  \riuch  allows  specifications  to  be  developed  which  will  result  in  known  fire  behavior  in  full-scale 
applications. 

5.4  Generalization  of  the  Current  Upward  Flame  Spread  Model 

The  computer  model  developed  to  predict  upward  flame  spread  and  burning  of  a  wall  has  been 
designed  to  be  generalized  to  allow  modeling  a  fire  against  a  wall  as  well  as  to  predict  comer  fire  scenarios. 
This  requires  the  implementation  of  additional  heat  transfer  modeling  involving  the  new  fire  source  as  well 
as  wall  to  wall  heat  transfer. 
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6.0  CONCLUSIONS 


A  computer  model  has  been  developed  that  successfully  addresses  upward  fire  growth  on  vertical 
surfaces.  The  agreement  with  experimental  results  is  good.  Assumptions  have  been  made  based  on  the 
limited  information  for  flame  heat  transfer  rates  to  cause  burning  and  spread,  and  more  complete  experimental 
results  are  needed.  Small  variations  in  the  material  properties,  flame  height  correlation,  and  flame  heat  flux 
after  ignition  cm  have  a  si^ficant  effect  on  fire  growth.  Variations  in  ignition  and  flame  spread  properties 
affect  the  predictions  significantly  as  well.  Therefore,  material  properties  must  be  accurately  determined  to 
predict  the  flame  spread  characteristics  of  the  material.  Despite  the  sensitivity  of  the  model  to  many  input 
parameters,  the  model  demonstrated  good  results  over  a  range  of  experimental  conditions. 

The  predicted  heat  release  rates  for  PMMA  in  both  the  small  and  full-scale  wall  fires  are  in  good 
agreement  with  the  data  using  two  different  flame  height  correlations.  The  use  of  two  flame  height 
correlations  in  the  model  is  due  to  the  absence  of  side  walls  in  the  experiments  conducted  in  reference  (1 1), 
heat  flux  distribution,  and/or  width  effects.  ’ 

Approximate  wall  heat  flux  profiles  have  been  developed  based  on  the  steady-state  experimental  data 
available  in  the  literature  (refer  Figure  3).  There  is  considerable  scatter  and/or  tmeertainty  in  the  data  above 
the  pyrolysis  zone  for  burning  wall  material.  The  fluxes  in  the  flame  region  have  been  measured  in  the  range 

of  20-40  kWnr  and  fall  rapidly  at  y  >  yj.  Additional  work  is  required  to  improve  correlations  of  these  heat 
fluxes. 


The  properties  of  the  materials  used  in  the  model  are  either  directly  found  finm  the  literature  or 
deduced  from  bench-scale  experimental  data.  Further  work  in  the  modeling  and  measurement  of  solid  pha^^ 
response  of  materials  to  heat  is  required  to  improve  prediction  capability. 

ft  is  the  intent  of  this  woric  to  present  a  simple  flame  spread  model  that  integrates  many  of  the  features 
and  capabilities  of  those  described,  ^t  requires  a  minimum  of  input  data  and  which  accurately  predicts  the 
fire  growth  along  vertical  walls  subject  to  ignition  sources.  Currently,  it  is  generalized  with  respect  to  the 
wall  material.  Ultimately,  the  model  will  be  generalized  so  that  a  number  of  two  and  three  dimensional 
geometries  may  be  modeled. 

The  ultimate  goal  of  these  modeling  efforts  to  predict  and  evaluate  the  fire  performance  of  U.  S.  Navy 
materials,  products,  and  assemblies  and  to  provide  a  link  between  bench-scale  ignition  and  heat  release  rate 
measur^ents  in  the  cone  calorimeter  with  full-scale  fire  performance  to  ensure  an  environment  safe  from 
destructive  fires  for  the  Navy  requirements.  To  achieve  this  goal,  this  model  will  be  used  to  predict  fire 
perfoimance  of  the  Navy  materials,  and  allow  specifications  to  be  developed  which  will  result  in  known  fire 
behavior  in  full-scale  applications. 

The  woiic  has  developed  a  correlational  method  which  relates  bench-scale  fire  test  results  to  full-scale 
comer  fire  test  results.  Fu^er  development  of  the  model  wiU  broaden  our  capabilities  to  correlate  bench- 
and  full-scale  results.  Ultimately,  a  full  predictive  link  between  bench-  and  full-scale  can  be  developed 
through  development  of  this  model. 
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Appendix  A 
Input  Data 


Integer  1-7;  Vegion  control  flag  (ciment  veision  reqdres  a  7  to  nm 

The  vereiOT  number  parameter  specifies  tiie  version  of  flame  sraead  model  for  which  tiie  inniit 
prepared.  This  prevent  diflferent  input  formate  to  be  used.  ^ 

fat^erO  or  1;  Burner  control  flag  0  -  normal  i  ■  neversM  off  ' 

^toegerOorl;  Debug on'off  Q-off  i.pn  ^ — 

Debug  on  printe  various  values  and  diagnosed  inessages  to  several  files  and  send  additional  date  to  the  scre^ 


iiteger  1-4:  Analysis  type 

^^eral  analj^is  flags  available  in  the  model.  There  can  be  a  line  fire  source  across  the  entire  lower 

i  1  of  a  wall,  and  or  pool  fire  c4>sure  4^ 

Lto  the  wall.  I- Lme  burner,  2  -  Partial  Ime  burner,  3  -  Pool  fire,  and  4  -  Pbol  fi!^  com^ 

titeggr  0  or  1;  Cone  calorimeter  test  1  -  tine  0  - 

Cone  test  exposure  has  no  ignition  soutceand  has  a  unifonn  flux  over  the  surtece  area. 

Real  >  0;  Height  Height  ^ifies  the  height  of  the  vertical  sample  in  meto  " 

jteal>0;  Wdflg  Width  specifies  the  width  of  the  vertical  sample  in  meters 
Real>0;  Tliiclaiess;  Tludmess  specifies  the  thickness  ofthesarnple  in  inetets 

Ihtegen  Nunilierof  nodes  in  tire  horizontal  diii»rti>yn 
liitegen  NunAerof  nodes  in  the  veitiad  diiection 

spr^  on  the  vertical  sample  by  breaking  up  the  surtece  into  a  Ian® 
number  of  elements  and  each  element  is  bounded  by  four  ntvti^  up  a  targe 

Real;  Ihennal  iigitia  (hpc)  (kWirf  ly  sec 

The  thermophysical  property  of  the  vertical  surface  fiom  open  literatere. _ 

Real:  Ignition  ten^iatuie,  K 

Ignition  temperature  of  tiie  material  fi-om  opailiterahire 

Real:  Heat  of  comhBtion  (ki1<^ 

Heat  of  combustion  of  die  material  fiom  open  literature. 


A-1 


Real  0-1:  Fractional  of 


consumed 


mass 

Specifies  a  fiaction  of  mass/area  that  is  involved  in  the  burning  process.  This  may  represent  char  or  non- 
combustible  constituents. 


Real  >  0:  Density, 

Density  of  the  material  fiom  open  literature. 


Real:  Ityrolysis  tempeiatme,  K 
Available  in  the  open  literature. 


Real:  Ambient  tempeiatme,  K 

Ambient  temperature  is  the  temperature  of  the  ambient  atmosphere. 


Real:  Initial  ten^ratme,  K 

Typically,  the  initial  temperature  is  the  same  as  ambient  temperature,  but  in  some  cases,  for  the  ignition 
simulation  die  sample  is  preheated  at  preset  tenqierature  fixim  external  radiation  source,  then  die  initial 
temperature  is  different  ^m  ambient  temperature. 


Integer  1-4:  Fire  type  (flame  height  coireiation) 

Flame  height  correlation  for  a  line  fire  against  a  wall  or  wall  fire  described  in  Section  22. 


Red>0:  Inclosed  heat  flux  (kW/mP) 

External  heat  flux  fiom  a  radiant  panel  or  cone  heater. 


Real  >  0:  Flame  heat  flux  after  ignition  (kW/nf) 

Experimentally  determined  incident  heat  flux  fiom  wall  flames  to  the  sur&ce  desoibed  in  the  Section  22. 


Real  >  0:  Tenfietature  difTerence,  K 

Specifies  temperature  rise  of  material  above  initial  tempoature  where  the  surfiice  temperature  calcnlarinns 
switches  fiom  Scant  method  to  Runge-Kutta  method.  Usually  0.5  to  1.0  K. 


Ihl^er  1-4:  Flem  heat  flux  correlation 

Upon  the  ignition  of  the  wall,  the  luminosity  of  the  flame  may  increase  ovct  a  typical  of  mediane  ignitor 
flame.  To  allow  for  this,  the  model  allows  ^e  use  of  different  heat  flux  coireiation  after  ignition  of  the  wall 
as  described  in  Section  22. 


Real  >  0;  Beak  heat  flig,  used  only  with  flame  heat  #  4 

The  experimental  peak  heat  flux  from  flame  of  burning  material. 

IntegerO  -  2;  0  -  do  not  use  special  line  fire  fluxes,  1, 2  -  me  high  and  low  respectively. 
1  or  2  specified  that  a  material  heat  flux  exposure  be  used  for  the  line  burner. 

Real;  Convection  factor  (W/mF-IQ 

The  convection  factor  used  when  calculating  heat  loss  fiom  cone  test  analysis. 


Real  0 -  1;  Emissivity/Ahsoiptivi^ oftematerial 

Radiative  emissivi^/Absorptivity  of  material,  used  for  computing  heat  loss  without  gain. 

Real>l;  IMs  is  the  irimnam  ratio  between  the  fire  height  and  the 

This  presents  the  flame  height  from  falling  below  the  pyrolysis  height  for  low  heat  release  rate  materials. 


Real  >  0:  Time  max,  sec 


The  time  is  die  lengdi  of  time  over  which  the  simiilatifm  takes  place. 


Real  >  0  <  Time  max;  Calculation  inteival,  sec 


^i^on^°"  interval  is  the  time  interval  between  eadi  writing  of  tiie  ou^jut  to  the  final  results  of  die 


Real>0;  HRR/unit  leigtli,  kW/m 
Line  burner  strength. 


Red  >  0;  Time  to  bumerto  rise,  sec 

Allows  line  burner  to  rise  fixwn  0  to  final  value  over  a  finite  time. 


Integer^O;  Numberof voids 

Any  air  gaps  or  inert  material  on  a  vertical  surface  left,  represented  by  voids. 
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A  - 1  -  Input  data:  Small-scale  PMMA 

7  Version  control  flag  (current  version  requires  a  1  to  run) 

0  Burner  control  flag  0  -  normal  1  -  never  shut  off 

0  Dd)ug  on 

1  Analysis  type 

0  Cone  Calorimeter  test  1  -  true  0  •  false 

0.90  Height  (m) 

0.20  Width  (m) 

0.0254  Thickness  (m) 

2  Number  of  nodes  in  the  horizontal  direction 

500  Number  of  nodes  in  the  vertical  direction 

0.60  Thermal  inertia  (kpc)  (kW/m^  K)^  sec 

593.0  Ignition  temperature  (Tjg  )  (K) 

25000.0  Heat  of  combustion  (kJ/kg) 

1.00  Percent  of  the  mass  burned  of  (1-  this  is  what  remains) 

1200.0  Density  (kg/m^) 

650.0  Pyrolysis  temperature  (K) 

298.0  Ambient  temperature  (K) 

298.0  Initial  Temperature  (K) 

2  Fire  type  (Hame  height  correlation) 

0.0  Inqx)sed  flux  (kW/m^) 

21.0  Flame  flux  aftCT  ignition  (kW/m^) 

2.0  Ten^erature  difference  (°C  or  K) 

2  Flame  heat  flux  correlation 

0.0  Used  only  with  flame  heat  flux  correlation  #  4 

0  0  -  do  not  use  special  line  fire  fluxes,  1, 2  -  use  forms 

0  Convection  factor  (W/m^-K) 

0.9  Bnissivity/Absorptivity  of  the  material  surface 

1.1  This  is  the  minimum  ratio  between  the  fire  height  and  the  pyrolysis  height 

1500  Time  max  (sec) 

1.0  Calculation  interval  (sec) 

1.0  HRR/unit  length  (kW/m) 

60.0  Time  to  burner  to  rise  (sec) 

0  Number  of  voids 
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A- 2  Input  data:  Full-scale  PMMA 

7  Version  control  flag  (current  version  requires  a  1  to  run) 

0  Burner  control  flag  0  -  normal  1  -  never  shut  off 

0  Debug  on 

1  Analysis  ^pe 

0  Cone  Calorimeter  test  1-true  0- false 

5.0  Height  (m) 

0.58  Width  (m) 

0.0254  Thickness  (m) 

2  Number  of  nodes  in  the  horizontal  direction 

500  Number  of  nodes  in  the  vertical  direction 

0.60  Thermal  inertia  (kpc)  (kW/m^  K)^  sec 

593.0  Ignition  temperature  (Tjg )  (K) 

25000.0  Heat  of  combustion  (kJ/kg) 

1 .00  Percent  of  the  mass  burned  of  ( 1  -  this  is  what  remains) 

1200.0  Density  (kg/m^) 

650.0  Pyrolysis  temperature  (K) 

298.0  Ambient  temperature  (K) 

298.0  Initial  Temperature  (K) 

3  Fire  type  (Hame  height  correlation) 

0.0  Imposed  flux  (kW/m^) 

3 1.0  Flame  flux  afto  ignition  (kW/m^) 

2.0  Temperature  difference  (®C  or  K) 

2  Flame  heat  flux  correlation 

0.0  Used  only  with  flame  heat  flux  correlation  #  4 

0  0  -  do  not  use  special  line  fire  fluxes,  1, 2  -  use  forms 

0  Convection  factor  (W/m*-K) 

0.9  Emissivity/Absorpdvityofthe  material  surface 

1.1  This  is  the  minimum  ratio  between  the  fire  height  and  the  pyrolysis  height 

1500  Time  max  (sec) 

1.0  Calculation  interval  (sec) 

1.0  HRR/unit  length  (kW/m) 

60.0  Time  to  burner  to  rise  (sec) 

0  Number  of  voids 
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Input  data:  Full-scale  Plywood  (Test  #  1) 

Version  control  flag  (current  version  requires  a  1  to  run) 

Burner  control  flag  0  -  normal  1  -  never  shut  off 

0  E)ebug  on 

1  Analysis  type 

0  Cone  Calorimeter  test  1  -  true  0  -  false 

2.4  Height  (m) 

0.61  Width  (m) 

0.0127  Thickness  (m) 

2  Number  of  nodes  in  the  horizontal  direction 

500  Number  of  nodes  in  the  vertical  direction 

0.475  Thermal  inertia  (kpc)  (kW/m^  K)^  sec 

623.0  Ignition  temperature  (Tjg )  (K) 

1 5000. 0  Heat  of  combustion  (kJ/kg) 

0.64  Percent  of  the  mass  burned  of  ( 1  -  this  is  what  remains) 

473.0  Density  (kg/m^) 

700.0  Pyrolysis  temperature  (K) 

305.0  Ambient  temperature  (K) 

433.0  Initial  Temperature  (K) 

2  Eire  type  (Flame  height  correlation) 

5.2  Impost  flux  (kW/m^) 

37.0  Flame  flux  after  ignition  (kW/m^) 

1 .0  Temperature  difference  (°C  or  K) 

6  Flame  heat  flux  correlation 

0.0  Used  only  with  flame  heat  flux  correlation  #  4 

0  0  -  do  not  use  special  line  fire  fluxes,  1,2-  use  forms 

0  Convection  factor  (W/m^-K) 

0.5  Emissivity/Absorptivity  of  the  material  surface 
1 .08  This  is  the  minimum  ratio  between  the  fire  height  and  the  pyrolysis  height 

300.0  Time  max  (sec) 

1 .0  Calculation  interval  (sec) 

0.25  HRR/unit  length  (kW/m) 

15.0  Time  to  burner  to  rise  (sec) 

0  Number  of  voids 
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A- 4  Input  data:  Full-scale  Wood  Particle  Board 

7  Version  control  flag  (current  version  requires  a  1  to  run) 

0  Burner  control  flag  0  -  normal  1  -  never  shut  off 

0  Debug  on 

1  Analysis  type 

0  Cone  Calorimeter  test  1  -  true  0  -  false 

1.8  Height  (m) 

0.30  Width  (m) 

0.013  Thickness  (m) 


2  Number  of  nodes  in  the  horizontal  direction 
500  Number  of  nodes  in  the  vertical  direction 
0.475  Thermal  inertia  (kpc)  (kW/m^  K)^  sec 


623.0  Ignition  temperature  (Tig )  (K) 

14000.0  Heat  of  combustion  (kJ/kg) 

0.65  Percent  of  the  mass  burned  of  (1-  this  is  what  remains) 

600.0  Density  (kg/m^) 

700.0  Pyrolysis  temperature  (K) 

298.0  Ambient  temperature  (K) 

298.0  Initial  Temperature  (K) 

5  Fire  type  (IHame  height  correlation) 

0.0  Impost  flux  (kW/m^) 

30.0  Flame  flux  aftCT  ignition  (kW/m^) 

1.0  Temperature  difference  (°C  or  K) 

1  Flame  heat  flux  correlation 

0.0  Used  only  with  flame  heat  flux  correlation  #  4 

0  0  -  do  not  use  special  line  fire  fluxes,  1, 2  -  use  forms 


0 

0.9 

0.0 


Convection  factor  (W/m^-K) 
Ermssivity/Absorptivity  of  the  material  surface 


This  is  the  minimum  ratio  between  the  fire  height  and  the  pyrolysis  height 


3000.0  Time  nuix  (sec) 

1.0  Calculation  interval  (sec) 

22.0  HRR/unit  length  (kW/m) 

20.0  Time  to  burner  to  rise  (sec) 


0  Number  of  voids 
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